Chained structure of directed graphs with applications to social and transportation networks

The need to determine the structure of a graph arises in many applications. This paper studies directed graphs and defines the notions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell$$\end{document}ℓ-chained and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\ell ,k\}$$\end{document}{ℓ,k}-chained directed graphs. These notions reveal structural properties of directed graphs that shed light on how the nodes of the graph are connected. Applications include city planning, information transmission, and disease propagation. We also discuss the notion of in-center and out-center vertices of a directed graph, which are vertices at the center of the graph. Computed examples provide illustrations, among which is the investigation of a bus network for a city.

A fundamental topological property of a graph, which will be briefly recalled in the section "Notation and some properties of graphs and networks", is multipartivity. The nodes in an m-partite graph can be split into m disjoint subsets V i , i = 1, 2, . . . , ℓ , called partite sets, with connections occurring only between the subsets, but not within the subsets. When ℓ = 2 , the graph is said to be bipartite. A refinement of bipartivity for undirected graphs, referred to as the chained structure of the graph, was introduced in Concas et al. (2021). The chained structure characterizes undirected multipartite graphs; an m-chained graph has only edges between nodes that belong to "subsequent" partite sets V i and V i+1 , i = 1, 2, . . . , ℓ − 1 (and vice versa). This paper extends the notion of chained graphs from undirected graphs, to directed graphs. The chained structure reveals the "depth" of a graph, i.e., how many steps it may take to go from a specified node to any other node, by following edges along their direction.
In Concas et al. (2021), we used chained graphs to identify central nodes by introducing the position centrality measure for nodes of an undirected graphs. This notion is a generalization of closeness centrality. Central nodes are identified by their location in the chained structure. For an overview of other centrality measures; see Borgatti (2005), Estrada and Higham (2010), Estrada and Rodriguez-Velazquez (2005). This paper generalizes position centrality to directed graphs. Specifically, for directed graphs that have directed spanning trees, we define in-position and out-position centralities of a node by examining two different types of directed spanning trees associated with the graph; see Gabow and Myers (1978) for a discussion on directed spanning trees. These centrality concepts shed light on the ease of communication within a network. In the two sections devoted to numerical examples, we compare them to other existing centrality measures. In general, it is impossible to state which centrality measure is the best, as the concept of centrality takes different meanings in different applications. What we show is that position centrality, by varying the value of the parameter on which it depends, is able to spot specific aspects of a network that are not detected by traditional measures, and that depend upon the underlying chained structure.
The identification of the chained structure of a directed graph also can be useful for detecting the presence of anti-communities, i.e., node subsets that are loosely connected internally, but have many external connections with the rest of the graph. Several methods have been developed to identify anti-communities in undirected graphs; see Concas et al. (2020), Estrada and Knight (2015), Fasino and Tudisco (2017). The relation between clustering and community detection in directed graphs has been discussed in Laenen and Sun (2020). In Concas et al. (2021), we illustrated how the chained structure may be used for introducing a density measure for computing an "anti-community score" for undirected graphs. We extend this measure to directed graphs in the present paper. To the best of our knowledge, while the identification of anti-communities has been studied in the literature (Estrada and Knight 2015) (see also Concas et al. 2020;Fasino and Tudisco 2017), the identification of near-anti-communities (which are associated with a small anti-community score) has not been discussed yet. This paper is organized as follows. The section "Notation and some properties of graphs and networks" introduces notation and discusses general properties of graphs that will be used later. Directed chained graphs are defined in the section "Directed ℓ-chained graphs and their adjacency matrices". They can be studied with the aid of directed spanning trees. This is discussed in the section "Directed chained graphs and directed spanning trees". The chained structure naturally leads to the concept of position centrality, defined in "Position centrality and some applications" section. Nodes with the largest position centrality are referred to as central nodes. Some data sets deriving from real-world applications, including a social network, are analyzed in the section "Some examples". The section "A case study about position centrality" sheds light on the properties of center nodes by considering a case study concerning a bus transportation network. Finally, the section "Conclusion" contains concluding remarks.

Notation and some properties of graphs and networks
is a set of nodes, or vertices, and E = {e i } m i=1 a set of edges, which connect the nodes. Two nodes v i and v j , for i = j , are said to be adjacent if there is an edge from node v i to node v j . In this context, an undirected edge between the nodes v i and v j points both from v i to v j and from v j to v i . The node v i is said to be connected to the node v j if there is a path from v i to v j , that is, if there is a sequence of edges {e r s } k s=1 such that e r 1 originates from v i , e r k points to v j , and if e r s points to v ℓ , then e r s+1 starts from the same node for s = 1, 2, . . . , k − 1 . A cycle is a path that starts and ends at the same node v i .
An undirected graph is connected if each pair of distinct nodes is connected by a path. A directed graph is said to be strongly connected if for each vertex pair (v i , v j ) the node v i is connected to the node v j , and the node v j is connected to v i . A directed graph is said to be semi-connected if for each vertex pair (v i , v j ) either the vertex v i is connected to the vertex v j , or v j is connected to v i . A directed graph is weakly connected if there is a path between each vertex pair (v i , v j ) in the underlying undirected graph, that is, in the undirected graph obtained by replacing all directed edges by undirected ones. We refer to Estrada (2012) and Newman (2010) for discussions on graphs and their properties.
An unweighted graph G with n vertices can be represented by an adjacency matrix A = [a ij ] n i,j=1 with a ij = 1 if there is an edge from vertex v i to vertex v j ; otherwise a ij = 0 . Since an undirected edge can be thought of as being made up of two directed edges (in opposite directions), the adjacency matrix of an undirected graph is symmetric; the adjacency matrix of a directed graph is nonsymmetric.
Multipartivity and, in particular, bipartivity are fundamental topological characteristics of graphs that model interactions between different types of objects. Bipartite graphs contain vertices that can be partitioned into two disjoint vertex subsets V 1 and V 2 , such that there are no connections between vertices in the same subset. Assume that the n vertices of a bipartite graph G are separated so that the first n 1 vertices make up the vertex set V 1 and the remaining n 2 = n − n 1 vertices make up the vertex set V 2 . Then the adjacency matrix A of G is of the form where C 1 ∈ R n 1 ×n 2 , C 2 ∈ R n 2 ×n 1 , and O denotes a zero-matrix of suitable order. If the graph G is undirected, then C 2 = C T 1 , where the superscript T denotes transposition. (2.1) In undirected ℓ-chained graphs, the nodes are divided into ℓ disjoint subsets so that there are edges only between nodes belonging to "adjacent" node sets, that is, all edges from a node in V i point to a node in V i+1 or in V i−1 for some i. This kind of partitioning is discussed in Concas et al. (2021).

Directed ℓ-chained graphs and their adjacency matrices
The directed chained graphs introduced in this section generalize the notion of undirected chained graphs defined in Concas et al. (2021).
2), such that v i ∈ V 1 and all edges from vertices in the set V j point to vertices in the set V j+1 for j = 1, 2, . . . , ℓ − 1 , where the chain length ℓ is the largest number of vertex subsets V j with this property. The vertex subset V j+1 is said to be adjacent to the vertex set V j .
The chain length ℓ of a directed ℓ-chained graph may depend on the choice of the initial vertex v i . After a suitable permutation of the nodes, the adjacency matrix A of a directed ℓ-chained graph G = {V, E} becomes upper block bidiagonal with zero diagonal blocks, where the submatrix A i ∈ R n i ×n i+1 describes the connections from vertices in V i to vertices in V i+1 , for i = 1, 2, . . . , ℓ − 1.
Example 3.1 Consider the graph of Fig. 1. This is a 3-chained graph with the chained node sets V 1 = {v 1 , v 2 } , V 2 = {v 3 } , and V 3 = {v 4 } . The initial node can be chosen to be either v 1 or v 2 . The adjacency matrix is where we can choose the submatrices Assume that a graph is known to be directed ℓ-chained for some ℓ ≥ 1 , but that the value of ℓ is not known. Moreover, let a permuted version of the matrix (3.1) be known (for some unknown value of ℓ ). Thus, the available adjacency matrix is of the form where P is a permutation matrix that modifies the vertex ordering. Given the adjacency matrix A , we are interested in determining the vertex subsets V 1 , V 2 , . . . , V ℓ in Definition 1, as well as the number of sets ℓ ≥ 1 . A method for determining if a directed graph is ℓ-chained and partitioning the nodes into subsets is described by Algorithm 1. Given an adjacency matrix A of a directed graph, the first node subset V 1 is obtained by considering the column indices j such that A ij = 0 for each row index i; see line 1 of the algorithm. Then the other vertex subsets are determined by identifying the blocks in A that describe connections with nodes in the preceding node subset (line 6). If it is not possible to determine the first vertex set, or if during the process it results that some node is connected to a vertex in a preceding subset, then the graph is not ℓ-chained. This process gives a constructive proof of the following result.
The definition of directed ℓ-chained graphs is quite restrictive. To be able to discuss properties of a larger set of directed graphs, we relax the requirements of Definition 1 to allow edges between vertices in the vertex subset V i to vertices in vertex subset V j for some j ≤ i with j not much smaller than i.

Definition 2
The directed graph G = {V, E} is said to be directed {ℓ, k i }-chained with initial vertex v i if it has the chained structure described in Definition 1 with the extension that edges from vertices in the set V j are allowed to point to vertices in the sets V max{j−k i ,1} , . . . , V j , V j+1 for j = 1, 2, . . . , ℓ − 1 and some k i ≥ 0 . The integer k i , which we refer to as the lower bandwidth, is the largest integer with this property.
We note that Definition 1 corresponds to the situation when k i = −1 for all i in Definition 2.

Definition 3
The minimal lower bandwidth, k, of a directed chained graph is defined as where the minimum is over all initial vertices v i in the vertex set V ⊂ V that gives maximal chain length ℓ . When k is the minimal lower bandwidth, the graph is said to be directed {ℓ, k}-chained.
The {ℓ, k}-chained structure is quite general. We conjecture that any weakly connected graph with n nodes is {ℓ, k}-chained for some n ≥ ℓ > k ≥ −1 . A small value of k indicates that information in the graph flows in a preferred direction, with small back propagation. This structure can be investigated by means of spanning trees as described in the section "Directed chained graphs and directed spanning trees".
Example 3.2 Consider the directed graph G shown in Fig. 2. It is a directed {5, 2} -chained graph with initial vertex v 1 . If one removes the edge from vertex v 4 to v 2 , the graph becomes a directed {5, 1}-chained graph with initial vertex v 1 . If one continues by removing the edge from v 3 to v 2 , then a directed 5-chained graph with the same initial vertex is obtained.
The adjacency matrix analogous to Eq. (3.1) for a directed {ℓ, k}-chained graph G = {V, E} can be represented by a lower block Hessenberg matrix when the nodes are suitably ordered. Here the block A ij represents edges that point from the vertex subset V i to the vertex subset V j . All superdiagonal blocks A i,i+1 are nonvanishing, because if all entries of the block A i,i+1 were zero, then there would be no edges from the vertex subset V i to vertices in the subset V i+1 . But this would contradict the fact that the graph G is directed {ℓ, k}-chained.
If the minimal lower bandwidth, defined by Eq. (3.2), is k = 0 , then there is at least one edge from a node to another node in the same vertex subset. The adjacency matrix corresponding to such a graph is upper block bidiagonal when the nodes are suitably ordered. Similarly, a lower bandwidth k = 1 indicates that when the nodes are suitably enumerated, the adjacency matrix can be represented by a block tridiagonal matrix. More generally, a small lower bandwidth (3.2) indicates that there only are edges between vertex subsets V j with close indices.
The following result shows that for strongly connected directed {ℓ, k}-chained graphs, directed cycles will be observed if k ≥ 1 . For semi-connected or weakly connected directed graphs, cycles are not guaranteed to exist.
Proposition 2 Let G = {V, E} be a strongly connected directed {ℓ, k}-chained graph with vertex partition V = V 1 ∪ · · · ∪ V ℓ . Assume there are no edges between vertices belonging to the same vertex set and that k ≥ 1 . Let e j,i ∈ E represent a directed edge from vertex v j to v i , where v i ∈ V i and v j ∈ V i+s for 1 ≤ s ≤ k . Then there exists at least one directed cycle that starts at v i , contains the edge e j,i , and ends at v i . The possible minimum length of the directed cycle is s + 1.
Proof Since the graph G is strongly connected and there are no edges between any nodes in the same vertex subset, the shortest possible directed path from vertex v i to v j has length s as shown below where v i t ∈ V i+t for t = 1, 2, . . . , s − 1 . Combining this path with the edge e j,i determines a directed cycle of length s + 1 .
Identification of the {ℓ, k}-chained structure of a directed graph (if present) sheds considerable light on properties of the graph, including the presence of anti-communities. Anti-communities are vertex subsets W i , i = 1, 2, . . . , q , of V such that there are many fewer edges from nodes in W i to nodes in W i , than from nodes in W i to nodes in W j for j = i . For instance, the node subsets V j of an ℓ-chained graph are anti-communities. Recent discussions on anti-community detection for undirected graphs can be found in Concas et al. (2020), Estrada and Knight (2015), Fasino and Tudisco (2017). There are several methods and measures that allow one to identify communities or clusters, such as the intra-cluster density which, for undirected graphs, is defined as the ratio of the number of internal edges and the number of all possible internal edges; see Fortunato (2010). An analogous density measure for computing the anti-community score for undirected graphs was introduced in Concas et al. (2021). Here, we extend this measure to directed {ℓ, k}-chained graphs.

Definition 4
The anti-community score ρ ∈ [0, 1] for a node subset V i of the node set V of a directed {ℓ, k}-chained graph is the ratio of the number of directed edges between the vertices in V i and the total possible number of directed edges between them. An anticommunity with score ρ is said to be a ρ-anti-community.
We remark that the anti-community score aims at identifying an approximate anticommunity as a node set for which ρ takes a small value. A large value of ρ does not necessarily identify a community, because it does not consider the connections between the nodes in V i and those not contained in V i .

Example 3.3
For directed ℓ-chained graphs with node subset partitioning (2.2), the subsets V i , for i = 1, 2, . . . , ℓ , are 0-anti-communities, because there are no internal edges. For a directed {ℓ, k}-chained graph described in Definition 2, the subset V i has a positive anti-community score ρ i when it has internal edges. If ρ i is small, then the subset V i may be considered as an approximate anti-community.

Directed chained graphs and directed spanning trees
The chained structure of a spanning tree T for an undirected graph G is used in Concas et al. (2021) to determine a chained structure for a graph G , if such a structure exists, and to approximate a graph without a chained structure by a graph with such a structure. In this section, we consider directed graphs that have directed spanning trees. We remark that not all directed graphs have a directed spanning tree. The directed spanning trees are employed to partition the node set V into subsets V i that determine directed ℓ -chained graphs; cf. (2.2). This approach to partition the node set V is applied to partitioning node sets of directed graphs that have a directed spanning tree, but do not possess a chained structure, and provides an approach to approximate a directed graph G without chained structure by a directed graph with chained structure.
We first briefly review results for undirected graphs. Let G = {V, E} be an undirected graph. A spanning tree for G is a subgraph T = {V, E ′ } that is a tree and contains all the vertices of G ; see, e.g., Concas et al. (2021), Deo (1974), Newman (2010). A spanning tree T is not uniquely determined by G and, in particular, depends on the chosen initial vertex of the tree, the so-called root.
When the graph G is directed, two different types of spanning directed trees, the outtree (or arborescence) and the in-tree, can be defined; see Deo (1974). We will employ both these directed trees.
Definition 5 An out-tree rooted at node v i for a directed graph G = {V, E} is a subgraph T i out = {V, E ′ } of G that is a tree with the same vertices as G , and such that for every vertex v j , for j = i , there is only one directed path starting at v i and ending at v j in the tree.
that is a tree with the same vertices as G , and such that for every vertex v j , for j = i , there is only one directed path from v j to v i in the tree.
In an out-tree, information may flow from the root to each vertex in the graph, while in an in-tree information may flow from any vertex to the root. In the first case, the root is a good source of information for the nodes of the graph; in the second case, the root is a good receiver.
Out-trees and in-trees exist for every vertex of a directed graph only if the graph is strongly connected. Any vertex in a semi-connected graph belongs to an out-tree or an in-tree. This follows from Proposition 3 below. We remark that this property is not guaranteed to hold for a weakly connected graph.
Proposition 3 Let G = {V, E} be a semi-connected directed graph. Then the graph G has at least one out-tree and one in-tree.
Proof Let v i , v j ∈ V be arbitrary distinct vertices. Then either v i is connected to v j , or v j is connected to v i . Assume there is a directed path P from v i to v j . If all the vertices of G except for v i and v j are on the path P, then P is an out-tree rooted at v i and an in-tree rooted at v j .
Let u be a vertex of G that is not on the path P. Assume that there is neither a directed path from u to v i nor a directed path from v j to u; otherwise, we extend P by including u as a root. Then E contains directed paths from v i to u and from u to v j . Therefore an outtree rooted at v i and an in-tree rooted at v j are obtained.
Proposition 4 Let G be a directed graph. If the vertex v i of G is the root of both an outtree T i out and an in-tree T i in of G , then the graph G is strongly connected.
Proof Let v i satisfy the assumption of the proposition. Then for any vertex v j , j = i , there is a directed path from v i to v j and vice-versa. Hence, for every pair of vertices (v k , v j ) , k, j = i , there is a directed path from v k to v j passing through v i and vice versa. It follows that the directed graph G is strongly connected.
Each directed spanning tree has a directed ℓ-chained structure (2.2). For out-trees, the root of the tree is the only vertex in the first set V 1 of the chained structure, and the partition of the vertex set V is determined by the relation between the vertices of the tree. Thus, the vertex set V 2 contains the children of the root and, in general, the vertex set V i contains the children of the vertices in V i−1 , i = 2, 3, . . . , ℓ . For an in-tree, the root belongs to the last vertex set V ℓ , and the partition of the vertex set is determined by following the direction of the edges backwards until one reaches the set V 1 , which contains the leaves farthest away from the root.
A directed out-tree for G does not exist, for example, if the first set V 1 contains more than one node. Similarly, a directed in-tree does not exist when the last set V ℓ contains more than one node.
The process of generating directed spanning trees for a directed graph is illustrated in the following example.
Example 4.1 Consider the directed graph G shown in Fig. 3. It is semi-connected. The out-tree T 1 out and the in-tree T 3 in rooted at v 1 and v 3 , respectively, are displayed in Fig. 4. These are the only out-trees and in-trees for the graph G . Their directed chained structure is illustrated in Fig. 5.
The chained structure of a directed spanning tree T of G can be used to detect, or approximate, the directed chained structure of G . The chained structure of G might not be unique, as it depends on the starting vertex and the directed spanning tree T .
Definition 7 Let T = {V, E ′ } be an out-tree (or in-tree) for the graph G . A directed ℓ -chained vertex set decomposition for T is said to be a directed ℓ-chained vertex set decomposition for G . We will refer to leaves of T as leaves of G. Now consider the in-tree T 3 in in Fig. 4. In this case, We have already mentioned that some semi-connected graphs may not allow an ℓ -chained partitioning for an arbitrarily chosen initial vertex. For example, vertex v 2 in Example 4.1 neither can be the root of an out-tree nor of an in-tree that spans the graph.
Let D = E \ E ′ be the set of the edges in G that are not in T , and let C(T ) denote the graph obtained by adding the edges in D to the spanning tree T . The graph C(T ) coincides with G and inherits the chained structure of T . The directed chained structure of the spanning trees T 1 out and T 3 in in Fig. 4 e j ∈ D there is an index 2 ≤ i ≤ ℓ − 1 such that e j connects a vertex in V i to a vertex in V i+1 .
If G is compatible with T , then the graph G = C(T ) is directed ℓ-chained. If instead there is at least one edge connecting a vertex in V i to a vertex in V i−k , for i = k + 1, k + 2, . . . , ℓ , and k ≥ 0 is the maximal number with this property, then the graph G = C(T ) is directed {ℓ, k}-chained.
The graphs C(T 1 out ) and C(T 3 in ) , obtained by adding the missing edges to the spanning trees of Fig. 5, are displayed in Fig. 6. The former graph is {4, 1}-chained and the latter one is {3, 1}-chained.

Position centrality and some applications
The notion of position centrality for vertices of an undirected network was introduced in Concas et al. (2021). It is a generalization of closeness centrality. This section generalizes position centrality to directed graphs by defining the in-position and out-position centralities of a node. The in-closeness centrality of a node measures how close this node is to those it is receiving information from, while the out-closeness centrality of a node shows how close the node is to the nodes it is sending information to.
Let (#V i ) denote the number of vertices in the set V i .
Definition 9 Let us assume that an out-tree T out = {V, E ′ } rooted at the node v for the directed graph G exist. Moreover, let V 1 , V 2 , . . . , V ℓ be the directed ℓ-chained structure, starting at vertex v, determined by the tree. For a fixed p ∈ R , the out-position centrality of v is defined as 3 Fig. 6 The directed graphs C(T 1 out ) and C(T 3 in ) corresponding to the directed spanning trees in Fig. 5. The edges in D = E \ E ′ added to the trees are drawn in red We refer to a vertex v c with the smallest out-position centrality as a p-out-center vertex.
Definition 10 Let us assume that an in-tree T in = {V, E ′ } rooted at the node v for the directed graph G exist. Moreover, let V 1 , V 2 , . . . , V ℓ be the directed ℓ-chained structure, ending at vertex v, determined by the tree. For a fixed p ∈ R , the in-position centrality of v is defined as We refer to a vertex v c with the smallest in-position centrality as a p-in-center vertex.
The in/out-position centralities depend on the spanning tree chosen. They can be defined for every node only if the directed graph is strongly connected. The out-center vertex can be described as an "information transfer station", such that it can "easily" send information to all the other vertices in the graph. A similar interpretation holds for the in-center vertex, which acts as an information sink. The following example illustrates how the in/out-position centralities of a vertex can be computed by using the chained structures starting from the vertex.
We turn to the in-position centrality of vertex v 2 . Consider the in-tree rooted at v 2 with vertex set partitioning V 1 = {v 3 , v 4 } , V 2 = {v 1 , v 5 } , and V 3 = {v 2 } . We have Since the graph G is strongly connected, we can compute the in/out-position centralities for all the other vertices similarly. When p = 1 2 and p = 1 , the vertex v 3 has the smallest out-position centrality. This indicates that v 3 is the out-center vertex. The in-center vertices are v 2 and v 5 for p = 1 2 and p = 1 . When p = 5 , the out-center vertices are v 1 and v 5 , while the in-center vertex is v 1 .
P out 1 (v 3 ) = 1 · 2 + 2 · 2 = 6, P in 1 (v 2 ) = 1 · 2 + 2 · 2 = 6, P in 1/2 (v 2 ) = 4.24, P in 5 (v 3 ) = 96. In a semi-connected directed graph, the vertices can be divided into three subsets: O , which contains vertices connected to every other vertex in the network, I , whose elements are vertices to which every vertex can send information, and M , which contains intermediate vertices. There may be a non-empty intersection between the sets O and I . Out-position centrality is defined only for vertices in O , while in-position centrality can be computed for vertices in I . Since every vertex belongs to at least one spanning tree, semi-connected graphs are directed ℓ-chained or directed {ℓ, k}-chained.
Vertices for weakly connected directed graphs also can be divided into the above three subsets O , I , and M . However, the sets O and I may both be empty, since out/in-trees are not guaranteed to exist. Hence, weakly connected directed graph may not possess a chained structure.

Some examples
This section describes a few examples concerned with directed graphs. For each graph, we analyze the presence of anti-communities by identifying its directed chained structure. The out/in-center vertices are identified by computing the smallest out/in-position centralities. Knowledge of the chained structure is beneficial in the following contexts: • Information dissemination in a social network: we are interested in determining directed ℓ-chained or directed {ℓ, k}-chained structures with initial vertex (center vertex) such that information from this node can reach all other individuals in the least amount of time, where we assume that the time is proportional to the path length. Similarly, we may be interested in determining which individual(s) can collect information from all other vertices in the least amount of time. Moreover, in a directed {ℓ, k}-chained graph, the presence of an edge e j,i from vertex v j ∈ V j to vertex v i ∈ V i , for i < j , indicates the possibility of feedback of the information from v j to v i . The minimal lower bandwidth k shows the minimal length of the path from v j to v i . • Prevention of the spread of an infectious disease: let the edges of a directed chained graph represent the spread of an infectious disease among subjects that are represented by nodes. An edge e j,i from vertex v j ∈ V j to vertex v i ∈ V i , for i < j , represents a secondary infection of v i from v j . It is reasonable to prevent the spread of disease by detecting and possibly eliminating the out-center vertex. In the context of COVID-19, it is important that out-center vertices be vaccinated. Similarly, it can be important to protect an in-center vertex from infection from other nodes. Vaccination may be one way to achieve this.
To graphically illustrate the chained structures revealed by the model discussed above, we first consider the following two small directed graphs: • ibm32 (32 vertices, 126 edges): collected from the IBM 1971 conference advertisement. After removing self-loops, the graph has 94 edges. It is available at https:// sparse. tamu. edu/ HB/ ibm32/.
These networks are not social networks, but nevertheless will be seen to have structure that can be studied with the concepts introduced in the present paper. The network n2c6b10 is represented by a weighted graph. We consider the corresponding unweighted graph obtained by setting all weights to 1. Figure 8 displays the out-tree T 20 out (left) and the in-tree T 20 in (right) for the ibm32 network. Both trees are rooted at vertex v 20 of the graph and have maximal chained structure length ℓ = 7 . The set C(T 20 out ) , which contains the additional edges that are not in T 20 out , is shown on the left-hand side of Fig. 9. It has {ℓ, k}-chained structure with minimal lower bandwidth k = 4 , that is, edges in C(T 20 out ) from vertices in the subset V i are allowed to point to vertices in the subsets V i−4 , . . . , V i , V i+1 for i = 5, . . . , ℓ . The righthand side of Fig. 9 displays the graph C(T 20 in ) with minimal lower bandwidth k = 4 . We conclude that the graph ibm32 is directed {7, 4}-chained with initial vertex v 20 . Consider the chained structure determined by the in-tree T 20 in . There are four 0-anti-communities (the first subset V 1 contains only one vertex) and 3 anti-communities with scores ρ 3 = 0.10 , ρ 4 = 0.11 , and ρ 5 = 0.06 . Moreover, since both the out-tree and in-tree are rooted at vertex v 20 , it follows from Theorem 4 that the graph is strongly connected. Starting from each vertex, an out-tree and an in-tree are constructed and their associated chained structures are determined. Figure 10 displays the chain length ℓ and the lower bandwidth k of the {ℓ, k}-chained structure associated with the out-tree (left) and the in-tree (right) rooted at each vertex of the graph ibm32. The property of "exact length" in the legend represents the maximal chain length and the property "not exact length" indicates that the chain length is not maximal. The symbol • in the figure displays the chain length of each structure as a function of the initial vertex v j , j = 1, 2, . . . , 32 . When the chain length is maximal, we use the symbol .
The symbols * and × display the lower bandwidth k j of the chained structure with initial node v j for j = 1, 2, . . . , 32 ; see Definition 2. The lower bandwidths of the chained structure with the maximal chain length are displayed by * symbols; for the other chained structures, the symbol is × . Among the lower bandwidths associated with the maximal chain length, the smallest k j is the minimal lower bandwidth. Hence, Fig. 10 shows the maximal chain length to be ℓ = 7 and the minimal lower bandwidth k = 4 . It can be seen that the maximal chain length and the minimal lower bandwidth are not achieved for each starting or ending vertex.
The left-hand side of Fig. 11 displays the 1-out-position centrality, i.e., the out-position centrality for p = 1 , of each vertex of graph ibm32. The 1-in-position centrality for each vertex is shown on the right-hand side of Fig. 11. The out-center vertices of the graph are v 2 and v 3 , and the in-center vertex is v 10 . Figure 12 displays the out-tree T 28 out with maximal chain length and its corresponding graph C(T 28 out ) for the graph n2c6b10. We note that this graph is directed {4, 0}-chained with initial vertex v 28 . Since the out-tree T 28 out is the only spanning tree of n2c6b10, the initial vertex is the out-center vertex and the graph is semi-connected. The graph n2c6b10 has three 0-anti-communities and the node subset V 2 is an anti-community with ρ 2 = 0.04.
We now determine the directed chained-like structure and center vertices for the medium-sized directed graph gre_1107 with 1107 vertices and 5664 edges. This graph arises from simulation studies in computer systems and is available at https:// math. nist. gov/ Matri xMark et/ data/ Harwe ll-Boeing/ greno ble/ gre_ 1107. html. We refer to the graph as gre. After removing self-loops, this graph has 4557 edges. The out-tree T 808 out and in-tree T 644 in with maximal chained structure length are displayed on the left and the right of Fig. 13, respectively. The graph C(T 808 out ) is {53, 4}-chained and the graph C(T 644 in ) has a {53, 5}-chained structure. The chain length and lower bandwidth of {ℓ, k}-chained structures starting and ending at each vertex of the graph gre are shown in Fig. 14. Only one directed out-tree, T 808 out , and one directed in-tree, T 644 in , are found to have maximal chain length. Their corresponding {ℓ, k}-chained structures have minimal lower bandwidth k = 4 and k = 5 , respectively; this is not visible in the figure because of the density of the symbols. Hence, the graph gre is a directed {53, 4} -chained graph with initial vertex v 808 . It has five 0-anti-communities and 48-anticommunities with the minimal score ρ = 0.01 and maximal score ρ = 0.5 . The graph gre is strongly connected since each vertex has both out-trees and in-trees.
The 1-out-position centrality and 1-in-position centrality of each vertex of the graph gre are shown in Fig. 15. The out-center vertex is v 400 with 1-out-position centrality P out 1 (v 400 ) = 10129 , and the in-center vertex is v 7 with 1-in-position centrality P in 1 (v 7 ) = 9030. For the previous test networks, of small to medium dimension, it was possible to determine both a spanning out-tree and an in-tree. Now, we investigate the presence of  • twitter (3656 vertices, 188,712 edges) available from http:// wiki. gephi. org/ index. php/ Datas ets, reproduces the connections of some part of the Twitter social network; • wikivote (8297 vertices, 103,690 edges). The nodes in the network represent Wikipedia users and a directed edge from node i to node j represents that user i voted on user j in an administrator election (Leskovec and Krevl 2014); • gnutella (10,879 vertices, 39,994 edges) is the p2pGnutella04 network from Leskovec and Krevl (2014); • foldoc (13,380 vertices, 120,700 edges) is an on-line searchable dictionary http:// foldoc. org/: an edge from term i to term j exists in the network if in the FOLDOC dictionary the term j is used to describe the meaning of term i. The network is available at http:// vlado. fmf. uni-lj. si/ pub/ netwo rks/ data/. • math (13,840 vertices, 195,330 edges) is available from Leskovec and Krevl (2014) and represents the interactions on the stack exchange web site Math Overflow https:// matho verow. net/. In particular, a direct edge is present between node i and node j if user i commented on user j's answer. Three of the above networks admit either a spanning out-tree or an in-tree. Two of them do not, so we eliminated the out/in dangling nodes, that is, vertices that do not have incoming edges or outgoing edges, respectively. This pre-processing is reflected in a different number of nodes in columns 3 and 5, than in column 2. The absence of dangling nodes is a necessary condition for the existence of an out/in spanning tree, but not sufficient, as the results in Table 1 confirm. The results show that there are real-world networks resulting from important applicative settings that have a directed spanning tree and that inherit from it a chained structure.
We now analyze in more detail the network twitter. To illustrate the different center nodes determined by varying the value of p in Definitions 9 and 10, we analyzed this graph for the p values reported in the first column of Table 2. It turns out that no vertex admits an in-tree, while most of the nodes (3485) have an out-tree. The table shows that different out-center vertices are identified when p varies, even if there is some stability for p between 0 and 1. Both the depth ℓ and the minimal lower bandwidth k of the corresponding {ℓ, k}-chained structure can be seen to grow with p.
In the same table, we also report the out-center nodes identified by other wellknown centrality measures. The out-degree, betweenness centrality (Newman 2010), PageRank (Page et al. 1999), and hubs score from HITS (Kleinberg 1999), have been computed by the centrality function of Matlab. The hub-centrality (Benzi et al. 2013) has been computed by the hubauth package, developed in Baglama et al. (2014) and available at https:// bugs. unica. it/ cana/ softw are/. Table 2 confirms the well-known fact that centrality measures often disagree, making it hard to judge which result is the best. The table also points out that trees rooted at the vertices with largest position centralities tend to identify chained structures with a smaller depth ℓ and bandwidth k than trees rooted at nodes considered important with respect to other measures. Table 2 illustrates that when using position centrality, we are able to identify chained structure with the smallest ℓ and k values. Figure 16 displays the spanning trees rooted at the first two out-center nodes of Table 2. It is evident that the tree corresponding to the larger value of p, shown on the right, produces a chained structure for which the cardinality of the V j sets with small j is larger at the beginning of the sequence than for the tree on the left; the cardinality of the sets in the tree on the left are more balanced. However, it is difficult to understand the effect of this parameter on the choice of the center nodes without knowledge of the identity and history of the people defining the vertices. The transportation network analyzed in the next section aims to clarify the meaning of p-outcenter nodes, as well as to compare position centrality to other centrality measures.
To gain some insight into how the centrality measures considered in this section are related, we computed the Kendall rank correlation coefficient between the in/out position centrality with p = 1 (IPC/OPC) and degree (IDC/ODC), closeness (ICC/ OCC), betwenness (BC), hub/authority (HC/AC), and PageRank (PRC) centralities. The comparison has been performed on the ibm32, n2c6b10, gre, twitter, and bus-ca networks. The last network will be discussed in the next section.
The results are displayed in Fig. 17. The networks n2c6b10 and twitter do not appear in the graph on the left, because they have no in-center nodes; there is only one outcenter node in n2c6b10, so for this data set the comparison is meaningless, as it is clear in the graph on the right. As it is expected, the Kendall coefficients of position centrality and closeness centrality for both incoming and outcoming connections are 1, meaning that the agreement of the two ranks is perfect. Position centrality is seen to be strongly correlated to degree and PageRank for some of the data sets, but the two graphs demonstrate that the considered centrality indexes represent different features of the networks.

A case study about position centrality
To investigate the effect of the parameter p on the choice of the center vertices determined by the out-position centrality P out p (v) and the in-position centrality P in p (v) of a vertex v, we studied a transportation network for which it is possible to judge by common sense the results of the analysis. Like all transportation networks, it is closely related to the social behaviour of the individuals living in the area of interest.
We considered the bus network that serves the metropolitan area around the town of Cagliari in Sardinia, Italy. The area is about 65 km 2 , hosts 4.2 · 10 5 people, and includes the town of Cagliari as well as four smaller municipalities very close to Cagliari, contiguous in some parts: Monserrato, Selargius, Quartucciu, and Quartu Sant'Elena; see Fig. 18.
The bus network was constructed using data available on the web. We refer to this network as bus-ca. There are 970 bus stops. They define nodes. The distance between the bus stops is not available. We therefore measure distance as the number of bus stops The graph on the left concerns the incoming connections, the one on the right the outgoing ones. The computation has been performed for all the networks considered in the paper between the starting and ending nodes on the shortest path. The bus routes define edges. The resulting network is unweighted. Some bus routes depend on the direction of travel, e.g., because some streets are one-way. The bus network therefore is directed. A bus network, like most geographical networks, is strongly influenced by the landscape and urbanization. The Cagliari commercial center is located on the south-west border of the network, in front of the harbor. A large pond, the Molentargius Saline Regional Nature Park, is located in the center of the urban area. It separates the four municipalities from Cagliari, and prevents straight travel between them.
We computed the p-out-center and the p-in-center vertices of the bus network for different values of p. The results are reported in Table 3, where the bus stops are identified by their name. The table also reports the center nodes according to the degree, betweenness (Newman 2010), PageRank (Page et al. 1999), HITS (Kleinberg 1999, and hub/ authority centrality (Baglama et al. 2014;Benzi et al. 2013). Since betweenness and Pag-eRank do not distinguish out-centers from in-centers, only one center node is reported for them. For each in/out-center vertex, the depth ℓ (i.e., the distance between the center vertex, which is the root, and a most distant leaf of the spanning tree) is reported. Here the distance is measured in terms of the number of edges on the shortest path between the root and the leaf. We also report the minimal lower bandwidth k of the {ℓ, k}-chained graph structure.
The center nodes corresponding to p = 0.5 and p = 1 are in the commercial center of Cagliari, the part of town where the largest number of shops and restaurants are located, and where a large number of bus routes converge. The out-tree rooted at the out-center corresponding to the "San Benedetto" bus stop is displayed in the right pane of Fig. 19. The density of nodes in the upper part of the tree shows that many of the first chained sets V i , i.e., sets with small index i, contain a large number on nodes. This indicates that it is possible to reach a large number of destinations within a small number of bus stops, i.e., in a small time.
When p is significantly smaller than one, the out-and in-centralities defined in Definitions 9 and 10 give a smaller weight to sets V i with large cardinality than when p is larger than one. The effect is that the spanning tree rooted at the corresponding out-center node is more balanced; see the picture on the left of Fig. 19. This corresponds to a less rapid decay in the number of elements of the sets V i and also to a smaller depth of the tree. The center node "Legnano" is located in the Pirri district of Cagliari, while "Zuddas", and "Riu Mortu" are in Monserrato, a neighboring municipality. Both these zones join Cagliari with the small towns in the west part of the area, so they are barycentric for the network. It is possible to reach the farthest parts of the network from them after a relatively small number of bus stops. Table 3 Network bus-ca (970 vertices): we report the p-out-center and the p-in-center vertices for different values of p, together with the in/out-centers identified by other centrality measures For each vertex we also report the depth ℓ of the associated spanning tree and the minimal lower bandwidth k of the corresponding {ℓ, k}-chained structure

Fig. 19
Network bus-ca (970 vertices): on the left, the out-tree spanning the network rooted at the node corresponding to the "Zuddas" bus stop, the out-center vertex for p = −10, −1, 0.1 ; on the right, the out-tree rooted at the "San Benedetto" bus stop, out-center node for p = 0.5, 1 For a value of p somewhat larger than 1, say p = 5 , the center nodes are found in densely populated parts of Cagliari. When the value of p becomes very large, the center nodes are suburban bus stops that are served by "strategic" routes that connects them rather easily to the rest of the network. The depth of the spanning trees, that is, the maximal length of the routes starting from the tree root, increases monotonously with p. It is remarkable that the minimal lower bandwidth is rather large. This is a consequence of that there are some bus routes going back towards the center node with only a single bus stop before reaching the center.
The other centrality measures, reported in the lower part of Table 3, with the exception of the betweenness centrality, produce center vertices located in the Cagliari harbor area, the touristic center. In any case, these measures produce results complying with the traditional idea of centrality, while varying the parameter p in the position centrality gives the possibility to consider different aspects of this transportation network. The betweenness central vertex is difficult to interpret, as it is located in the central part of the Quartu Sant'Elena town, which does not appear to identify the real center of the network.

Conclusion
It is important to be able to identify interesting structural properties of directed graphs, because they shed light on how the vertices are connected. This paper introduces the notion of directed chained graphs and illustrates how it helps us to understand the structure of directed graphs. Also, the related notions of in-central and out-central nodes are defined and illustrated. The latter notions are quite intuitive and examples illustrate that they are helpful for identifying important nodes that differ from nodes that are identified by several popular available centrality measures.